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We study the formation and growth of laser-induced periodic surface structures (LIPSSs) with 
the finite-difference time-domain (FDTD) method. We use a recently proposed inter-pulse feedback 
method to account for the evolution of the surface morphology between each laser pulse sent to the 
surface of the processed material. This method has been used with an ablation-like mechanism, 
by removing material exposed to a light intensity higher than a given threshold. We propose an 
inverse mechanism, an expansion-like mechanism, able to grow structures that the ablation-like 
process cannot. This allows us to introduce the notions of constructive and destructive feedback 
and explains a strong contradiction between the standard Sipe-Drude theory and the experimental 
observations, i. e. the formation on metals of structures usually linked to wide band gap dielectrics. 


I. INTRODUCTION 

The study of laser-induced periodic surface structures 
(LIPSSs) has started with the first experimental observa¬ 
tion by Birnbaum et al^ They are created by irradiating 
the surface of solid material (dielectrics, semi-conductors, 
metals)Pl with an intense laser beam near the ablation 
threshold. LIPSSs are wavy nanometric structures char¬ 
acterized with both their orientation with respect to the 
incident wave polarization and periodicity A. A large va¬ 
riety of distinct LIPSSs morphologies^Hlol be formed 
depending on the material and the laser properties. 

We will consider two of these morphologies. The first 
type of structures to be considered are the most com¬ 
mon and also the first ones to have been observed. They 
have an orientation orthogonal to the incident laser beam 
polarization and a periodicity close to its wavelength, 
A ^ A. They are referred to as low spatial frequency 
LIPSSs (LSFLs) or as types because of the sinusoidal 
dependency between the angle of incidence and their pe¬ 
riodicity. 

The second type of structures have an orientation par¬ 
allel to the incident light polarization and a periodicity 
closer to A ^ A/Re(n), where n is the complex refractive 
index. Since the real part of the refractive index can vary 
significantly during intense laser irradiation, the period¬ 
icity of parallel structures can cover a wide range. They 
are sometimes referred to as LSFLs or high spatial fre¬ 
quency LIPSSs (HSFLs) depending on their periodicity, 
or type-d structures. 

The currently accepted theory is the Sipe-Drude 
theor}®Ki^ which accounts for the interaction between an 
incident plane wave and surface scattered waves caused 
by the material’s surface roughness. Sipe’s analytical so¬ 
lutions of Maxwell’s equations predict that a material 
with Re(n) > Im(n) (dielectrics) results in the formation 
of type-d structures and a material with Re(n) < Im(n) 
(metals or strongly ionized dielectrics) gives rise to type-s 
structures. 


Agreement is quite good with the experiments on sev¬ 
eral points. On wide band-gap dielectrics at low laser 
fluence, type-d structures are usually observed.^I^On nar¬ 
rower band-gaps or at higher fluence, the resulting struc¬ 
tures are rather of type-s because of the formation of 
dense plasma. On metals, type-s structures are observed 
most of the time, but in contrast to the Sipe-Drude the¬ 
ory, parallel structures have also been seen on met als an d 
this observation is currently without explanation 

Since the analytical method cannot predict the for¬ 
mation on metals of LIPSSs parallel to polarization, we 
turn to a numerical approach, the finite-difference time- 
domain (FDTD) methocP^®^ to investigate further. A 
feedback mechanism is also included to account for the 
growth of LIPSSs from one laser pulse to the next.^^^ An 
extension to the model allows us to add an important no¬ 
tion, the distinction between constructive and destructive 
feedback in the LIPSSs formation process. This turns 
out to be crucial for the formation of type-d structures 
on metals. 


II. NUMERICAL MODEL 

A. FDTD and discretization 

FDTD is a numerical method used to solve Maxwell’s 
equations and is based on their spatial and temporal dis¬ 
cretization. We first define a tridimensional domain, rep¬ 
resented in Figure in which half of the volume is occu¬ 
pied by the bulk of the material and the other half with 
void together with randomly distributed set of inhomo¬ 
geneities positioned at the interface. Each inhomogene¬ 
ity has the minimal volume allowed by the discretized 
space, i.e. one computational cell. Futhermore, 10% of 
the cells in the layer above the bulk are randomly filled 
with material. Each spatial boundary is completed with 
perfectly matched layers (PMLs), 20 cells wide, to avoid 
any non-physical boundary reflections. 
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Figure 2. Incident light is mostly directed to red spots. Am¬ 
plitude growth (constructive feedback) is achieved by deeper 
ablation under the surface minima. 


Figure 1. Geometry used in the FDTD calculations. Material 
is represented in orange, void in light gray and PMLs in black. 
PMLs are also present at vertical boundaries (not shown). 


All spatial dimensions are normalized with respect to 
the incident light wavelength and the spatial discretiza¬ 
tion is = A /20 in a domain of 10 A in x and y and 

2A in 2 ;, excluding PMLs. The time dimension is normal¬ 
ized with respect to the time of a complete optical period 
T and the temporal discretization is 5t = T/40, smaller 
than the maximal time step St^rnax allowed by the stabil¬ 
ity condition of the FDTD method {St^max = (20v^) 

in our case). A plane wave is then propagated from 
the highest layer of void towards the -^z direction and 
through the material roughness. Simulations run for 7 
full optical cycles and we extract the average intensity 
I{x^y^z) = {\E{x^y^ z)\'^) over the 3 final optical cycles. 


B. Drude model 


The optical properties of the material are calculated 
by the use of the Drude model to account for the dense 
plasma formation during the laser interaction.^^ The re¬ 
fractive index is defined as n = and the complex 
permittivity e is 
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where cj is the laser frequency, 7 is the collision frequency 
of the free carriers defined as the inverse of the Drude 
damping time td and Uy is the plasma frequency de¬ 
fined as ojy = e^Nejeorri^y^rrie. Finally, e is the electric 
charge, is the free carrier density, eo is the free space 
permittivity, m*p^me is the effective optical mass of the 
electrons and e is the material permittivity with no free 
carriers fixed at e = 4.84 in this work. 
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Figure 3. Incident light is mostly directed to red spots. Am¬ 
plitude growth (constructive feedback) is achieved by expan¬ 
sion of the surface maxima. 


C. Feedback 

To account for the growth of LIPSSs between two laser 
pulses, an inter-pulse feedback mechanism recently pro¬ 
posed by Skolski et al is included.^^^ After the first sim¬ 
ulation, the time averaged light intensity /(x, y, z) in the 
material is calculated and an isosurface s(x, y^ z) is traced 
along a threshold value of average intensity Ith‘ All the 
material above s{x,y,z) = Ith has received more inten¬ 
sity than the threshold and is thus removed before start¬ 
ing a new simulation with the processed surface. This 
ablation-like process is repeated as many times as needed. 
The threshold is chosen to remove a precise average thick¬ 
ness of material As after the first laser pulse and is kept 
constant for the subsequent pulses. As is also in units of 
A. 

Several distinct and complex mechanisms are responsi¬ 
ble for the modification of the surface morphology in ad¬ 
dition to ablation. For instance, photo-expansion will act 
differently on the surface because it will lead to expansion 
of the material where the light intensity is larger. We 
therefore introduce a second mechanism to account for 
any expansion-like mechanism with the same threshold 
method except that the surface modification is applied 
with the opposite sign. For example, if the difference 
between the actual surface position and the isosurface s 
at a certain point is 5 cells, an ablation-like process will 
remove 5 cells of material and an expansion-like process 
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Figure 4. Growth of periodic nanostructures with the inter-pulse feedback method. Subfigures (a) to (d) use param¬ 
eters (c^p/cj, 7 /ct;, As) = (1.7,1/16,-0.11) for pulse numbers 1,2,3 and 4 respectively. Subfigures (e) to (h) use pa¬ 
rameters (co^p/ct;, 7 /ct;, As) = (6,1,0.16) for pulse numbers 1,3,5 and 9 respectively. Subfigures (i) to (1) use parameters 
(cjp/o;, 7 /cj, As) = ( 6 ,1, —0.16) for pulse numbers 1,3,5 and 9 respectively. 


will add 5 cells. We use negative values of As to indicate 
an expansion-like mechanism. 

Several scenarios can occur between two inter-pulse 
feedback cycles. When LIPSSs are already formed on 
a surface, the calculated mean intensity I{x^y^z) can 
be maximal either at the surface minima (valleys) or at 
the surface maxima (peaks). If the energy is mostly de¬ 
posited at the valleys, two scenarios illustrated in Figure 

are possible. If the chosen surface processing mech¬ 
anism is ablation, the valleys get deeper and the total 
amplitude of the LIPSSs grows, leading to constructive 
feedback. In contrast, if the process is expansion-like, 
the valleys expand and the LIPSSs flatten as a result of 
destructive feedback. 

In distinction, if the energy is mostly deposited at the 
peaks, as illustrated in Figure the opposite effects are 
expected. The ablation mechanism will erase the LIPSSs 
in destructive feedback and the expansion mechanism 


will lead to amplitude growth as a result of constructive 
feedback. 


III. RESULTS & DISCUSSION 

To test these assumptions, we first apply 9 full inter¬ 
pulse feedback cycles on a material with dielectric op¬ 
tical properties (Re(n) > Im(n)) under the expansion¬ 
like mechanism. The parameters are (cj^/u;, y/cj. As) = 
(1.7,1/16, —0.11), corresponding to n = 1.4 + 0.064i and 
results are shown in Figure [^a)-(d). The analytical so¬ 
lution of the Sipe-Drude theory predicts the formation 
of structures parallel to the polarization, in agreement 
with this simulation. Furthermore, our calculation con¬ 
firms the iiiustration of Figure [^whereby structures grow 
from one iaser puise to the next when expansion is ap- 
piied, whereas they do not under abiation (not shown in 
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The next two sets of calculations consider materi¬ 
als with metallic properties (Re(h) < Im(h)). We 
again apply 9 full inter-pulse feedback cycles, but this 
time under the ablation-like mechanism with parame¬ 
ters (cjp/cj, y/cj, As) = (6,1, 0.16), corresponding to h = 
2.14 +4.21i. The simulation results are shown in Figures 
[^e)-(h) and are in agreement with the Sipe-Drude the¬ 
ory both predicting type-s structures. The constructive 
feedback under ablation mechanism is again qualitatively 
displayed in Figure 

Now, with the same three parameters, {ujpjuj, 

= (6,1,—0.16), we apply an expansion-like mechanism 
to obtain Figures |^i)-(l). One clearly sees that struc¬ 
tures with orientation parallel to the laser polarization 
are grown, a behavior not accounted for by the Sipe- 
Drude theory on metallic materials. This is in sharp 
contrast with the experimental fact that these types of 
structures are indeed observed on metals.f^^IIIl Qur sim¬ 
ulations suggest that the growth of type-d structures on 
materials with metallic optical properties is the result 
of energy deposition at the peaks with the ensuing con¬ 
structive feedback for the developement of the maxima 


as illustrated in Figure 

IV. CONCLUSION 

We have extended a recently proposed feedback-by¬ 
ablation methocP^ to account for expansion-like mecha¬ 
nisms in the LIPSSs formation process. This allows us 
to simulate the growth of a larger class of nanostructure 
morphologies, including parallel ripples on dielectrics and 
metals. 

We have thereby investigated the importance of con¬ 
structive and destructive feedback. In particular, we have 
proposed an explanation for the growth of parallel struc¬ 
tures (type-d) in the optically metallic regime, something 
not accounted for by the standard theory without feed¬ 
back. 
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